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Using variational trial wave function based on the string picture we study the motion of a single 
mobile hole in the stripe phase of the doped antiferromagnet. The holes within the stripes are taken 
to be = static, the undoped antiferromagnetic domains in between the hole stripes are assumed to 
have alternating staggered magnetization, as is suggested = by neutron scattering experiments. The 
system is described by the t — t' — t" — J model with realistic parameters and we compute the single 
particle spectral density. 
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I. INTRODUCTION 

The experimental discovery of static stripe-like struc- 
tures in some cuprate materials and nickelatesp has re- 
cently attracted considerable attention. This is particu- 
larly remarkable for the cuprate superconductors because 
stripe-like hole arrangements have been predicted k>ng 
ago by Hartree-Fock calculations for the^d-p modclEI as 
well as for the single-band Hubbard modeE - these would 
then in some sense be the first successful theoretical pre- 
dictions ever made for these materials. So far, there ap- 
pears to be experimental evidence for stripes in a CuC>2 
based material only for the concentration of x = 3D 1/8, 
where the material does not superconduct. This might 
suggest that the stripes in this case are a kind of CDW- 
instability, stabilized predominantly by Coulomb inter- 
action between holes, which occurs only at this special 
concentration. On the other hand, the rather precise 
scaling of the vector of magnetic incommensurability in 
\j&2-x Sr x Cu04 with x for hole concentrations x < 15% 
can be quite naturally explained by assuming that fluc- 
tuating stripes with an x-independent hole density of 1/2 
within the stripes and a mean spacing of ~ l/2x between 
stripes do exist, which are separated by Neel ordered do- 
mains with alternating staggered magnetization. 
In that sense, stripes could be a generic feature of a 
doped Heisenberg antiferromagnet. Exact diagonaliza- 
tion (ED) studies of the t-J model have shown evidence 
for stripe-like hole correlations in the ground stateo, as 
well as low-energy excited states where the holes form 
stripe-like arrangements!! These results cannot actually 
be compared to the experimental situation, because the 
relatively small size of clusters available for the ED tech- 
nique does not allow to accommodate a sufficiently large 
supercell as is observed experimentally. This is possi- 
ble, however, in density matrix renormalization group 
calculations and this technique also has produced strong 
evidenced for stripe-like hole arrangements in the ground 
state, or at least in very low excited states of the t-J 
model. 

In the present paper we do not want to address the ques- 



tion whether stripes really do form in the t-J model or 
what is the energy responsible for their stabilization. In- 
stead, we take the existence of a stripe-like hole arrange- 
ment for granted, and study the dispersion and spectral 
function of an additional hole, which is moving through 
the stripe 'background'. In other words: we want to dis- 
cuss the photoemission spectrum of a striped phase. This 
seems an interesting question, because recently angle re- 
solved photoemission (ARPES) data for La 2 - X Sr^CuC^ 
have become available!! This compound is known to 
show stripe order for x = 3-D 1/8 and the ARPES spectra 
show a rather nontrivial evolution from the hole disper- 
sion in the insulator (x — 3D0) which is very similar 
to Sr2Cu02Cl2P, via the formation of 'new states' near 
(•7T, 0) or (0, 7r) to a very different spectrum for the super- 
conducting phase (x ~ 0.15). Our intention is to check 
whether this can be reproduced by studying, in much 
the same way as the hole motion in an antiferromagnet 
is being studied, the motion of a single mobile hole in 
a static background of spins and immobile holes, shown 
in Figure |l|, which form a stripe pattern rather than the 
conventional undoped Neel state. 
We used the t-J Hamiltonian: 



H = 3D- 



H.c 



Si ■ Sj 



(1) 



Here 



At _ 



3£>ct (1- 



are the standard conditional 



electron creation operators, which avoid double occupan- 
cies, Si the operators of electron spin, and fii occupa- 
tion numbers. The hopping integral t{j between nearest 
neighbors was taken to be t = 3D1, that for 2 nd nearest 
neighbors was taken as t' = 3D — 0.3, the one between 
3 rd nearest neighbors was t" = 3D0.3. All other hopping 
integrals, as well as all exchange constants for other than 
nearest neighbors were taken to be zero. The nearest 
neighbor exchange constant was J = 3DQA. 
The hole motion in an antiferromagnet can be under- 
stood at least qualitatively within the string pictureBEj. 
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A hole created in a Neel state cannot propagate, because 
it leaves behind a trace of misaligned spins. This means 
that the magnetic frustration in the system increases to 
good approximation linearly with the distance travelled 
by the hole, whence the hole experiences an 'effective po- 
tential' which traps it around the site where it has been 
created. The wave function of the respective self-trapped 
state can be obtained by diagonalizing the Hamiltonian 
in the subspace of 'string states' where the hole is con- 
nected to its starting point i by strings of defects with 
various length and topology. True derealization then 
is possibly only by virtue of the transverse part of the 
Heisenberg exchange which can 'heal' the defects created 
by the hopping hole. In fact, if the first two defects cre- 
ated by the hole are repaired by the transverse exchange, 
this amounts to shifting the starting point of the string 
to a 2 nd or 3 rd nearest neighbor of i. The resulting next- 
nearest neighbor hopping band of width oc J is indeed 
well established by a variety of techniques, both numer- 
ically and analytically^. One may expect quite similar 
processes to occur also in a stripe-like background such 
as the one shown in Figure [j]. The holes in the stripes 
thereby are to be considered as static. Also in this case 
an added mobile hole cannot propagate freely, because 
it creates frustration in the remaining antiferromagnetic 
volume of the system. Again the spin-flip parts will be 
necessary to enable propagation across the = antiferro- 
magnetic stripes. The main difference then is that there 
may be additional processes when the hole passes through 
the hole stripes, and the fact that these hole stripes form 
something like 'A/2-plates' for the propagating hole, be- 
cause they separate antiferromagnetic domains with op- 
posite staggered magnetization. As already mentioned 
the latter property is an experimental constraint, which 
is necessary to explain the incommensurate low energy 
peaks in neutron scattering. The calculation for the 
stripe background thus will be somewhat more compli- 
cated than for the case of the pure Neel state, but in 
principle the formalism is precisely the same. In the fol- 
lowing we outline the respective calculation and present 
the results. As will be seen below, the computed single- 
particle spectra show some resemblance with the experi- 
mental results of Ino et at. 



II. SPIN POLARONS IN THE STRIPE 
STRUCTURE 

Our goal is to discuss the motion of holes in the stripe 
structure shown in Figure [lj In the following, we denote 
this state as |$o)- It corresponds to a concentration of 
static holes of 1/8= 3D12.5% and has been suggested by 
experimental The 32 sites contained in the rectangle in 
Fig.|l| form the unit cell of this structure. Some theoreti- 
cal studies suggest different structures!] but also indicate 
a tendency towards formation of stripes. Charge or spin 
order within the stripes has not actually been observed. 



The pattern depicted in Fig.|l| has been chosen for sim- 
plicity. The charge and spin fluctuations within stripes 
might, in principle, be taken into account in subsequent 
steps of approximation, but here we restrict ourselves to 
the simplest case possible, namely static holes and spins. 
We stress again that our aim is not to discuss the stability 
of the stripe structure, but rather to take it for granted, 
and concentrate on the properties of an additional hole 
injected into the stripe structure. We assume that its 
dynamics is represented by the standard t — J model, 
augmented by additional next-nearest neighbor (t') and 
next-nexttjaearest (t") hopping terms. Such terms are 
necessaryE3~t3 to explain the dispersion in the insulating 
Cu02 plane as measured by Wells et alB. 

As already mentioned, a mobile hole which has been 
created in an ideal antiferromagnetic background behaves 
as if it is trapped in a kind of potential wellEI. The reason 
is that its movement creates strings of misaligned spins 
in the antiferromagnetic structure and thus increases the 
magnetic energy. These strings of defects, tend to restrict 
the motion of the hole to the neighborhood of the cite 
where it has been created. The notion of this effective 
potential well for the hole is equivalent to the idea of 'spin 
bags' which was introduced in the context of the Hubbard 
modeO. Coherent motion of holes is possible due to 
relaxation of strings which is enabled predominantly by 
the transverse part of the Heisenberg exchange contained 
in the t — J model. The motion of a hole, then, has a 
twofold character. Because of t > J,t,t',t", the hole 
will oscillate rapidly (with a hopping rate oc in the 
neighborhood of the site to which it is attached by the 
string and only rarely (at a rate tx J -1 ) the exchange part 
will shorten the string of defects by two. This means that 
the starting point of the string, and hence the center of 
gravity of the oscillatory motion, is shifted by two lattice 
spacings. 

In order to describe this twofold dynamics of holes 
we first construct wave functions which describe a self- 
trapped state of the hole. We first introduce the opera- 
tors 

t <m> = 3£) ( c li c id + c h c hi) ■ ( 2 ) 

In an arbitrary spin background they shift holes by one 
lattice spacing. By consecutively applying the T/ it j\ to 
a the state c;|$o), we generate a family of string states. 
The strings of defects connect the hole to the site I. We 
denote such a string state by \l,Vi), whereby Vi stands 
for a set of numbers which give the topology of the string. 
We will choose V = 3D(ei, ea, . . .), where the unit vec- 
tor ej £ {±e x ,±e y } gives the direction of the i th hop 
which the hole has taken. For the wave function of a 
hole trapped in the vicinity of the site I we then make 
the ansatz: 

\* l ) = 3Dj2<*l,'P l \h'Pi). (3) 
{■Pi} 
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For the Hamiltonian we choose the sum of the kinetic 
energy and Ising part of the Heisenberg exchange: 



H = 3D-tJ2J2( l 



H.c. 



(i,3) 



(4) 



The kinetic energy has matrix elements of —t between 
two states \l,Vi) and \l,V[) if the respective paths can 
be transformed into each other by a single hop. The di- 
agonal matrix elements of H Q originate from the term 
JJ2{ij) {Si,zSj, z — riiUj/A). If we choose the state |$o) 
as the zero of energy for these diagonal matrix elements, 
the energy of any given string state is simply the number 
of additional 'broken' bonds times J/2. For example the 
energy of the state C(o,iVf |<&o) is ^/2, while for the state 
C(-i,i),tl$o) it is 3J/2. 

Once we truncate the subspace of string states (for exam- 
ple by taking into account only paths up to a maximum 
length) we can easily set up the matrix corresponding to 
Hq and diagonalize it numerically to obtain the ground 
state 



(5) 



This gives us the coefficients a/.-p, in (|3|) as well as the 
(/-dependent) ground state energy Ei. This procedure 
can be performed for any 'environment' which the hole 
has been created in. The variation of the wave function 
\^l) with I describes how the spin bag is 'deformed' by 
its local environment. This 'prediagonalization' of Hq 
reduces the number of degrees of freedom considerably, 
and thus makes the problem tractable even with the rela- 
tively large unit-cell for the stripe structure. This is also 
the main|-diffcrence compared to the related method of 
Trugmant2l, which obtains the ground state wave func- 
tion by constructing one Bloch states from each of the 
string states \l,V). Neglecting the excited states of Hq 
in the further steps of the calculation is justified provided 
the separation in energy between these excited states and 
the ground state is sufficiently large so that transitions 
into these states can be neglected. We have verified that 
the separation in energy between the ground state and 
the first excited state is in any case > J/2, in most cases 
even larger than J. Since the bandwidth for a single hole 
is « 2J, neglecting these excited states is probably not 
really an excellent but reasonable approximation. 
In our calculation we have used the part of the Hilbert 
space defined by the condition S z = 31? — 1/2, where S z 
is the z component of the total spin. The other pos- 
sible choice S z = 3D1/2 would give identical results, 
corresponding to a hole moving on the opposite sublat- 
tice. The condition S z = 3D — 1/2 may be fulfilled by 
removing an up spin from the stripe structure. There 
are 14 non-equivalent spin-up sites in the elementary cell 
in Figure The presence of the ferromagnetic bonds 
across the stripes brings about an extra complication: 



some of states that represent holes created in stripes are 
directly coupled by the 'fast' nearest neighbor hopping 
term. For example, the state 6(5 x),f |$o) may be obtained 

\by acting only once with the operator on cuu^\<f>o). 

/This cannot happen in the Neel state, where sites of like 
spin are always separated by at least 2 lattice spacings, 
and thus are not coupled by one nearest neighbor hop. 
To deal with this situation, we incorporate both states, 
C(5,i),f|^o) and c^iw^o), into one and the same self- 
trapped state \^i), which consequently is elongated in 
x-direction. In this way the number of non-equivalent 
starting points I for holes is reduced to the following 12 
sites: (0,1), (1,0), (1,2), (2,1), (2,3), (3,0), (3,2), (4,1), 
(5,3), (6,0), (6,2) and (7,3). Finally, it turns out that 
some starting points I for holes are equivalent in that the 
respective spin polaron states |\&;) may be transformed 
into one other by applying a symmetry transformation 
of the stripe structure. For example, the sites (0, 1) and 
(4, 1) which belong to the same elementary cell are equiv- 
alent. 

By means of computer algebra we can now create the 
states The algorithm thereby takes into account 

exactly strings of length up to 3 lattice spacings. For 
longer paths (which are less relevant) the approximation 
is made that their amplitude in I*;) does not depend on 
details of further steps. At this stage of the calculation 
we also make the assumption that the states \l,Vi) for 
different I and/or different Vi are always orthogonal. For 
longer paths this may not be true. The stripe pattern 
provides much more possibilities to reach identical states 
from different original configurations, not only by means 
of retracing paths as it is basically in the case of the Neel 
state, but we neglect this here. As already mentioned, 
any charge and spin fluctuations of the holes and spins 
which form the background stripe pattern are also ne- 
glected. 

It then turns out that the differences between the 
eigenenergies Ei of spin-polaron states for different ini- 
tial positions I of the hole are ~ J. The states which are 
lowest in energy thereby correspond to polarons centered 
on (0, 1), (4, 1) and equivalent positions. We see that de- 
spite the fact that the freedom of motion of a hole created 
at the site (0, 1) is strongly reduced by the (static) holes 
at (0,0) and (0,2), which inevitably results in a loss of 
kinetic energy, there is a net gain in energy compared to 
other sites. This is obviously due to the smaller num- 
ber of broken bonds created in the first few hops. This 
means first of all a lower exchange energy, and second a 
gain in kinetic energy because the hole can move more 
freely along the remaining paths. We do not pursue this 
any further, but we note that this may be one reason for 
the stability of stripes. 
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III. PROPAGATION OF SPIN POLARONS IN 
THE STRIPE STRUCTURE 

So far we have constructed wave function for self- 
trapped states at various positions in the stripe struc- 
ture. We now assume that the wave function for a hole 
added to the stripe structure is a coherent combination 
of these polaron states \^i). All processes which were 
neglected when constructing the self-trapped states |^;) 
will now be incorporated into an effective Hamiltonian 
H e ff which couples the \&i) centered on different I. As 
already mentioned, in doing so we implicitly assume that 
the excited states which come out of (||) are separated 
from the ground state by an energy which significantly 
exceeds the matrix elements of H e f f . 
The diagonal element (^fi\H e ff\^i) is essentially the 
eigenenergy Ei of the polaron-state It influences 

the probability that a particular site is occupied by the 
polaron. By means of computer algebra we have also set 
up additional diagonal and off-diagonal contributions to 
H e ff. Thereby we have confined the search to paths no 
longer than 2 lattice spacings. There are several hun- 
dreds of process which couple pairs of sites one of which 
(at least) belongs to the elementary cell. They form 29 
different categories of processes. We shall now discuss 
some examples. 

To begin with, the possibility of direct hopping to sec- 
ond nearest neighbors by virtue of the t'-term in the t — J 
model has not been taken into account when setting up 
the different Hq. Therefore, H e ff must take it into ac- 
count. Consider Fig.||. Shifting a hole created at the 
site (0,1) by (1,1) gives rise to a state that represents 
a hole created at (1,2). This process therefore couples 
the states |^(oi)) and l^n^))) the numerical value of 
the corresponding matrix element (^(0,1) l-f e/ / 1^(1.2)) is 
t'o!(Q i)Q!(i,2)' Here pairs of numbers in subscripts refer 
to position of holes, and the a's are obtained by diago- 
nalizing the Hq matrices. 

Next, Fig.^(a,b,c) depicts a process which contributes to 
the diagonal element for the spin-polaron state at (2, 1). 
After a single hop in the (1, 0) direction a hole created at 
(2, 1) occupies the site (3, 1). The amplitude of this state 
in |^(2,i)) is ai,T>i — 3-Da(2,i),(i,o)> where, as explained 
above, the first pair of numbers denotes the original po- 
sition of the hole whereas the second pair gives the di- 
rection of the hop. Next, let us assume that the t' term 
now shifts the hole to the site (2,2). This state has the 
coefficient Q!(2,i),(o,i) m 1^(2, 1)) whence we have found a 
contribution of i'ct(2,i).(i,o) a (2,i),(o.i) to the matrix ele- 
ment (* (2)1) |/r e// |*j 2 ,i.))' 

There are a few more types of processes that involve the t 1 
term. In Fig.||(d) we start with the string state obtained 
by creating a hole at (2, 3) and performing two hops in 
directions (1,0) and (0, 1). Shifting the hole in (—1, — 1) 
direction by virtue of t' , this is transformed into a state 
obtained by creating a hole at (3, 4) and shifting it in the 
directions (0,-1) and (-1,0) (Fig.|(e)). This type of 



process couples polarons at the sites (2, 3) and (3, 4) and 
the corresponding contribution to (^ r ( 2 ,3)|-ffe//|^ r (3,4)) is 
i'a(2,3),(i,o)(o,i) a (3,o),(o,-i)(-i,o)- Yet another process is 
possible due to the existence of parallel spins on near- 
est neighbor sites in the stripe structure. Shifting a hole 
created at the site (3,2) (see Fig.||(f)) by (1, 1) produces 
the same state as obtained = by two hops of this hole in 
the directions (0,1) and (1,0). The corresponding con- 
tribution to the matrix element {^(%^)\H e ff\^Ci,2}) is 

*' a (3,2)Q ! (3,2),(0a)(l,0)- 

The third nearest neighbor hopping term t" model gives 
rise to similar additions to H e f f . For example a contri- 
bution t"a( 5j3 )a( 7i 3) to (*(5,3)I-He//|*(7,3)> emerges due 
to the possibility of direct hopping = of a hole from 
(5, 3) — > (7, 3). Fig.|| shows a process whereby a hole cre- 
ated at the site (4, 5) (Fig.|](a)) which hopped three times 
in the directions (1,0), (0,-1) and (0, -1) (Fig.|(b)) 
is shifted by (0,2). The resulting state (Fig|3(c)) then 
represents a hole created at (5, 3) which hopped twice 
in direction (0, 1). That gives rise to a contribution 

- i ""(4,5),(i,o)(o,-i)(o,-i)"(5,3),(o,i)(o,i) to the matrix el- 
ement (\E r (4 5j I -f^e/ / 1 ^(5.3) ) ■ The minus sign stems from 
Fermi statistics and the definition of the operator Tuj\ 
which creates string states and the fact that the length 
of the first string is odd, while the length of the second 
is even. 

This concludes our discussion of the effects of t' and 
t" and we now proceed to a discussion of contributions 
to H e ff which originate from the exchange part of the 
t — J model. Figjj(a,b,c,d) depicts a process that in- 
volves pairs of parallel spins on nearest neighbor sites in 
the stripe structure. Fig.[|(a,b) shows holes which have 
been created at the sites (4,1) and (8,1). After three 
hops in the direction (—1, 0) the hole from the site (8, 1) 
reaches the site (5, 1). Fig.|](c) shows the corresponding 
string state which is a component of the spin polaron at 
the site (8, 1). By acting with the spin-flip term on the 
bond (6, 1) — (7, 1) this state is transformed into the state 
in Fig.^(d) , which in turn is a string state which is a com- 
ponent of the spin polaron at the site (4, 1). To be more 
precise, the state in Fig.^(d) can be obtained by a single 
hop in the (1, 0) direction from the state in Fig.^a. The 
matrix element describing the process in Fig.||(a,b,c,d) is 

(J/2)a( 4) i) ) (i ) o)a(8,i) ) (-i,o)(-i,o)(-i,o)- This mechanism 
of coupling spin polarons at the sites (4, 1) and (8, 1) also 
works for all longer string states where the hole in Fig.|] 
(c) and (d) has moved in the vertical direction. This mo- 
tion does not interfere with the action of the exchange 
term on the sites (6, 1) and (7, 1) which gives rise to the 
equivalence of string states formed by further motion of 
the hole in Fig.|](c) and (d). The above contribution to 
H e f f therefore should be supplemented by terms related 
to all relevant longer paths. 

Quite generally it is obvious that Figures ||,||, and [| rep- 
resent only a small fraction of all processes which may 
occur during the motion of a hole through the stripe 
background. There are various incquivalcnt positions of 
polarons within the unit cell and the different local ar- 
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rangement of spins and holes in the stripes gives rise to 
many non-equivalent processes which however may be 
easily found by means of computer algebra. 

As we have mentioned before, spin and charge fluc- 
tuations have been neglected in our calculation. They 
influence to some extent the = properties of the system. 
Nevertheless, it seems plausible as a first attempt to treat 
the stripe structure as a static object and discuss the = 
propagation of an excess hole in the stripe background. 
We expect to gain in this way some insight into electronic 
properties of stripes. 

In order to construct a coherent combination of spin 
polaron states we have to label the positions of polarons 
in a slightly different way. = We identify a site I by 
the vector R m which gives the position of the lower left 
corner of the elementary cell containing the site and the 
position n of the site within the elementary cell: 

l*i> = l*R m ,n) (6) 

A Bloch state |^) with momentum P then can be written 
as 

|i)=3i)-^e iR »> p ^|% Bin ) ) (7) 

where N' is a normalization constant, i.e. the number 
of elementary cells in a system with periodic boundary 
conditions. The momentum P has to be chosen from 
the reduced Brillouin zone [— f-] <£> [— f , f ]■ The states 
l^n) are normalized by definition: 

(* n \* n ,)=3D8 n , n , (8) 

The states |\£ n ') and |^„) for different sites n and n! in 
the elementary cell are orthogonal because up to the ac- 
curacy level of our calculation different string states are 
orthogonal. The matrix elements of the Hamiltonian for 
the coherent combinations \*5f n ) — then may be repre- 
sented in terms of the matrix elements for spin polarons 
which we already know: 

(*„|ff e// |* n ,) = 3^^ e iR m ' p (vl/ ,„| J ff e// |* Rm „„,)- 

m' 

(9) 

In order to get a feeling for the evolution with doping 
we have analyzed an analogous stripe structure, shown 
in Figure ^, which would correspond to the case of 1 
hole per 12 copper sites. The existence of stripes at this 
level of doping is as yet a speculation. Stripes, if they 
exist at all, may in this case have a dynamical charac- 
ter. Nevertheless, we have decided to use a structure 
which is analogous to the pattern that we used for the 
1/8 case. This can also give some indication as to how 
robust the various features in the calculated spectral den- 
sity are against possible changes of the structure. Finally, 
to conclude the discussion of the propagation of the mo- 
bile hole, we note that the dispersion relation for a hole 



in the originaLNeel state, can be carried out in precisely 
the same wayliS. We have also performed this calcula- 
tion to compare the dispersion for the Neel state and the 
stripe phase. 

IV. SPECTRAL FUNCTION 

We now proceed to discuss the one-particle spectral 
function and consider the term related to creation of an 
additional hole in the stripe structure. It can be written 
as 

A(k,c) = K*S +lfc) |cMl*o>| 2 « (« - - 4 S) )) , 

(10) 

where l^" 1 " 1 ^) is an eigenstate corresponding to an ad- 
ditional hole inserted into the stripe structure and in our 
calculation is identified with the solution of the eigenvec- 
tor problem for the Hamiltonian <^). The Ck,| operator 
may be written as 

Ck , T =3Z^]T e -* kR "=c„, T . (11) 

* n 

The summation here runs over all sites in the lattice and 
N = 3D32N' is = their total number. The v-th excited 
state is a linear combination of the states |^ n ): 

|4+ lfc ))=3£>53^ n) l*n>. (12) 

n 

The summation label n thereby refers to the 12 inequiva- 
lent sites in the elementary cell enumerated above, which 
may be occupied by a spin-up polaron. After some 
straightforward algebra = the relevant correlation func- 
tion takes the simple form 

($(+ 1,l >|c k , T |4 s) ) = 3D5 ModiPx+kx , 2n/8) , x 
*Wo« l (P, + * s ,W4),oyyE(^ T,) )* e ~ kR, <*''.»l c '.Tl*o>. (13) 

n,l 

The matrix elements (^/o.nlQ.tl = ^o) give the overlap 
between the bare hole created in the stripe structure = at 
site I and the self-trapped state centered on site n. Since 
the electron annihilation operator cannot create any spin 
defects, it can only produce the 'string of length zero'. 
Hence only the state c n ^\^o), which has the prefactor a n 
in I'I'o.n), can possibly have a nonvanishing overlap with 
Q iT |$ )- Therefore we have (*o,nKtl$o) = 3D5i, n a n . 
The only exception are the 'elongated' states along the 
ferromagnetic bonds. As was explained above, the states 
obtained by creating a hole either at (4, 1) or at (5, 1) 
are both included into the self-trapped state centered on 
(4, 1). Hence 



5 



(*(4,i),oI c J,tI $ o) = SDSm^a^i) - 5|,(5,i)Q!(4,i),(i,o)) 

(14) 

where the relative minus between the two contributions 
is due to Fermi statistics. 

We have employed the formula (|ll|), together with the 
above expressions for the matrix elements, in the eval- 
uation of the spectral functions. It should be noted 
that the resulting 'spectral weight' basically is a measure 
as to how well the respective many-body wave function 
matches the wave function of a photoelectron, which we 
model by a plane wave with momentum k. The actual 
spectral weight observed in an angle resolved photoemis- 
sion experiment may still differ from this: the spectral 
weight of a quasiparticle in the t-J model is known to 
have a significant k dependence due to the quantum spin 
fluctuations of the Heisenberg antifer-eamagneillj, in the 
more realistic singlec3 or three-bandEil Hubbard model 
an additional k dependence is brought about by charge 
fluctuations and/or interference between photoholes cre- 
ated on Cu and O sitesEil. Unfortunately these effects are 
beyond the scope of our approximation. 
On the other hand, the Brillouin zone of the supercell in 
Figure |l| is [—§,§] <8> [—f, f]- An angle resolved photoe- 
mission experiment with a momentum transfer of, for ex- 
ample, fc = 3D(tt, 0) thus actually probes the eigenvalue 
spectrum at (f,0). Any 'dispersion' which is measured 
on scales in k space larger than the extent of the reduced 
Brillouin zone therefore is in reality a 'spectral weight 
dispersion' and cannot be interpreted as a band disper- 
sion. This 'dispersion' then is essentially determined by 
the expression above, namely how well a photoelectron 
with the respective momentum matches the wave func- 
tion of the respective state within the 32-site unit cell. 
In that sense the above 'spectral weight' does very well 
have some significance. Therefore, one may expect that 
although the magnitude of the spectral weight obtained 
from (fi~3) cannot be directly compared to experiments, 
it still may give us a rough feeling as to where in (k, uj) 
space we may expect nonvanishing weight. 
The smallness of the reduced Brillouin zone brings about 
another complication: suppose we compute (as is usu- 
ally done) the spectral function at an equally spaced k- 
grid along some high-symmetry line, e.g. (0, 0) — > (tt, 0). 
Then, there are two possibilities: the step size of the k- 
grid may be commensurate with the size of the reduced 
Brillouin zone or not. To illustrate the first case, let us 
assume the step size to be tt/8. Then, we would actually 
be probing just two momenta in the reduced zone, namely 
(7r/8, 0) and (0,0). This would therefore give us an ar- 
tificial periodicity of the peak positions (although not of 
the spectral weight). On the other hand, if we choose an 
incommensurate step size (we choose tt jl as an example) 
we would walk through the reduced Brillouin zone like 
this: (0,0), (-f ,0), (f^), (-f ,0), (|f ,0)... that 
means more or less in a random sequence. Plotting the 
results for the spectral density in this way therefore is 
probably of very little significance. We have therefore 



decided to average the spectra for each momentum over 
a small neighborhood of the respective momentum and 
compute: 

A(k,u) = 3£>y^ w itj A h (k + iS x +jS y ,uj). (15) 

i,3 

We have chosen the weight function Wij to be a con- 
stant and used 5 x . y = 3D-g^e xy , and summed over 
—8 < i,j < 8. To some degree this also simulates the 
effect of the finite fc-resolution in an ARPES spectrum. 
Bearing this in mind we turn to the numerical results for 
the spectral density. Figures |6|-|To| then show the spectral 
function A(k, u>) calculated along various high-symmetry 
lines in the extended Brillouin zone. Distances in the 
k space between points in the sequence are equal. An 
important point for comparison with experiment is the 
following: we are using a hole language, which implies 
that the ground state of the hole is the state with the 
most negative energy. In a photoemission experiment 
the single-hole ground state then would actually be the 
first ionization state, i.e. the state with the lowest bind- 
ing energy, and would form the 'top of the band'. For 
comparison with the usual way of plotting experimental 
ARPES data, with binding energy increasing to the left, 
the energy axis in our figures would have to be inverted. 
To begin with, the direction (0,0) — > (tt, tt)) (see Figure 
fuj and (0,0) -► (0, 7r) (which has momentum parallel to 
the stripes, see Figure R) show a distribution of spectral 
weight which broadly follows the quasiparticle dispersion 
for the t — t' — t" — J model in the pure antiferromagnet 
(the latter is indicated by the full line). The direction 
(tt, 0) — > (tt, tt) (see Figure ||) might be interpreted as the 
superposition of two components, one of them following 
the antiferromagnetic quasiparticle dispersion, the other 
one forming a dispersionless band at « — 2.3t. It should 
be noted, however, that in this part of the extended Bril- 
louin zone the spectral weight of the quasiparticle band 
is practically zero in the antiferromagnet!, whence the 
dispersionless band is probably unobservable. The situa- 
tion is different for the direction (0, 0) — > (tt, 0) (i.e. the 
momentum is perpendicular to the stripes), see Figure]^. 
In this direction we see several more or less dispersionless 
'bands', and in particular a flat band of low energy states 
at rs — 2.3t. Since this range of momenta is within the an- 
tiferromagnetic Brillouin zone one may expect, based on 
the results for the antiferromagnet, that it has an appre- 
ciable spectral weight. A dispersionless band at w — 2.3i 
is also seen in the direction (0, tt) — > (tt, tt), but again thus 
is probably not observable because it is in the outer part 
of the extended Brillouin zone. Averaging the directions 
(0, tt) and (tt,0), which presumably corresponds to the 
experimental situation because the stripes do not have a 
uniform direction throughout the sample, one would thus 
get a 'broadened' version of the quasiparticle band in the 
insulator, plus a relatively flat band of new states, which 
actually form the first ionization states. It is quite tempt- 
ing to identify this with the experimental results of Ino et 
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at. At half- filling {x — 3D0), these data showed a disper- 
sion which is very much reminiscent of the 'quasiparticle 
band' observed in S^CuC^C^. As doping is increased, 
this band persists to some degree, although is getting 
more diffuse. In addition, around (tt, 0) a relatively dis- 
persionless band of 'new states' emerges, which forms 
the first ionization states. In order to clarify the charac- 
ter of the wave functions, we have also studied the real 
space wave function w n = 3D\j3v | 2 for different states 
(see Eq. (12)). This basically gives the probability to 
find the hole in the self-trapped state at site n within the 
unit cell. Figure |ll] shows this for a state which produces 
an intense peak at (it, 0), whereby the peak relatively far 
from the top of the band and thus would correspond to 
the 'remnant of the antiferromagnetic band'. It is ob- 
vious from the figure, that the mobile hole in this case 
avoids the hole stripes and resides mainly in the center 
of one of the antiferromagnetic domains. On the other 
hand, Figure [l^ shows the probability distribution for 
one of the 'new states' at (tt, 0). There the extra hole 
resides almost exclusively within the stripe. This shows 
again, that the new states, which appear upon doping, 
are related to states which are predominantly localized 
in the hole stripes. Of course one has to keep in mind 
that our approximation of static holes within the stripes 
will make the description of precisely these states the 
most questionable. Finally, Figure O shows w n for a 
state which produces a significant peak at (0, tt) which 
also belongs to the 'remnant of the antiferromagnetic 
band'. Again, we see that the mobile hole avoids the 
hole stripes, which is expected if the state somehow re- 
sembles properties of the undoped antiferromagnet. In 
order to understand the origin of some features observed 
in the spectra for the 1 /8 case we have also performed the 
calculation of the spectral function for the structure rep- 
resented by Fig. ^, which would correspond to the case 
of S = 3D1/12. Stripes at this level of doping, if exist, 
have a dynamical character and the applicability of the 
stiff structure represented by Fig. || is very questionable. 
Our aim was to check the robustness of some character- 
istics, rather than to derive any spectra comparable with 
experiments. Figures [l4| and [if] depict spectra calculated 
along the lines (0, 0) — > (tt, 0) and (0, 0) — > (tt, tt) respec- 
tively. We again notice a band at the energy « — 2.3t 
which emerges in the vicinity of (tt, 0) and disappears 
somewhere halfway to the point (0, 0). We have checked 
that this band may be attributed to states that repre- 
sent the extra hole which resides predominantly within 
the stripe. The rest of the spectrum seems to resemble 
the energy dispersion of a single hole in an antiferromag- 
net. We have not observed that for the 1/8 case when the 
size of antiferromagnetic domains is smaller. The states 
which form the 'antiferromagnetic band' correspond to 
the extra hole avoiding stripes. The spectrum for the di- 
rection (0,0) — > (i")7r)) however very broad, follows the 
quasiparticle dispersion of the pure antiferromagnet in a 
more evident way then for S = 3D1/8. These remarks in- 



dicate that the evolution of the spectral function between 
doping levels S = 3D0 and 5 = 3D1/8 may be attributed 
to formation of stripes. 

V. DISCUSSION 

In summary, we have calculated the dispersion for a 
single mobile hole in a static stripe structure or, put an- 
other way, a depleted Heisenberg plane, thereby using a 
reasonably realistic version of the t—t'—t"—J model. The 
strongest approximation presumably consisted in taking 
the holes in the stripes to be static, but we believe this 
is a reasonable first step in understanding the dynamics 
of holes in truly fluctuating stripes. In particular, the 
present calculation incorporates the effect of the oscil- 
lating direction of Neel order in the undoped antiferro- 
magnetic domains of the stripe pattern, as well as the 
influence of the important t' and t" terms on the hole 
motion. From the calculated dispersion and wave func- 
tion we have computed the 'spectral weight' as a function 
of energy. This spectral weight is basically the overlap of 
the wave function in the supercell with a plane wave of 
the respective momentum. In an ARPES experiment on 
a system with a relatively small reduced Brillouin zone as 
the stripe structure with its large unit-cell, it is basically 
this 'spectral weight dispersion' which is being measured. 
As expected, the primary influence of the stripes is the 
'broadening' of the spectral weight, which however still 
roughly follows the hole dispersion in the insulator, at 
least along some high-symmetry lines. This is to be ex- 
pected for a not too strong super-cell potential. In addi- 
tion, the stripes lead to the emerging of 'new states' near 
(tt, 0) at relatively low binding energy. More precisely, 
the spectral density in the stripe structure shows a rather 
sharp and dispersionless low energy band, which is pre- 
sumably most intense around (tt,0), and which actually 
forms the first ionization states. Analysis of the wave 
functions has shown, that the holes reside predominantly 
within the hole stripes in these new states, whereas the 
states responsible for the 'remnants of the antiferromag- 
netic band' in the spectral density have the holes pre- 
dominantly in the antiferromagnetic domains. It is in- 
teresting to note that this effect, namely the generation 
of low energy states at (tt, 0) by going from the undoped 
to the doped system, can also be observed in exact diaa- 
onalization calculations for the t — t' — t" — J modelE-X 
Whether this is due to formation of stripes even in the 
relatively small clusters studied remains to be clarified. 
In any case the effect of 'doping' (that means the modifi- 
cation of the Neel state to the stripe pattern with static 
holes and antiferromagnetic domains) on the spectral 
function in our calculation shows some similarityjvith re- 
cent ARPES measurements on Ij&2-x Sr x CuO40. They 
covered the range from an optimally doped supercon- 
ductor (x = 3-D0.15) to an antiferromagnetic insulator 
(x=3D0). It has been suggested on the basis of the in- 
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commensurate peaks in inelastic neutron scattering9@ 
that the family of LSCO systems reveals an instability 
towards a static spin-charge order of stripe form. The 
formation of these stripe structures is accompanied by a 
suppression of T c at a hole concentration 6 ~ 1/8. 
It should be noted that in the experimental ARPES spec- 
tra 'peaks' in the spectral density cannot usually be as- 
signed. The various 'bands' in the spectra become visible 
only upon forming the 2 nd derivative of the spectral den- 
sity with respect to energy. For the insulating compound 
La2Cu04 the 'dispersion' roughly resembles the known 
band structure of a hole in an antiferromagnet as it has 
been observed in S^CuC^C^: the top of the band ap- 
pears to be at (ir/2, tt/2), the band at (tt, 0) is at a signif- 
icantly (rj 200meU) higher binding energy. Upon doping 
this band structure of the insulator persists to some de- 
gree, but in addition 'new states' at (tt, 0) appear. These 
have a very low binding energy and indeed do form the 
first ionization states. Let us concentrate on the case of 
x = 3D0.12 which is closest to the configuration we have 
discussed. It is believed that stripe structures are perpen- 
dicular in nearest copper-oxygen planes. The spectra in 
the (0, 0) — > (0, tt) direction are a composition of spectra 
parallel and perpendicular to stripes. The measurements 
in that direction correspond rather well to a combination 
of our calculated spectra for the directions (0, 0) — > (0, 7f) 
and (0,0) — ► {tt, 0). The same holds true for the direc- 
tion (tt, 0) — » (tt, tt). Keeping that in mind we notice 
that the experimental and theoretical spectra show at 
least a qualitative agreement. In both cases we observe, 
along the directions (0,0) — > (tt, 0) and (0,0) — > (0, 7r), 
the 'splitting' of the spectrum and formation of new low 
energy states around (tt,0). This splitting is much less 
pronounced in the direction (0, tt) — > (tt,tt). The spectra 
along the direction (tt, 0) — > (tt, tt) are very broad. That 
prevents assignment of a band for some range of k points 
in that direction. Finally, the spectra along the direc- 
tion (0, 0) (tt, tt) are very broad and it is practically 
impossible to distinguish any band. 
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FIG. 1. Stripe structure for 12.5% of doping which have 
been applied = in the calculation of the hole dispersion. Solid 
circles depict immobile holes within the stripes, arrows repre- 
sent spins in antiphase AF domains that separate the stripes. 
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FIG. 2. Some processes related to direct hopping to next 
nearest neighbors. 

FIG. 3. A process related to direct hopping to third nearest 
neighbors. 



FIG. 4. A processes related to the exchange term. 
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FIG. 5. Structure applied in the calculation for the case of 
1 hole per 12 sites. 



FIG. 6. Spectral functions in the (0, 0) — » (tt, 0) direction 
for points separated by an equal distance. 



FIG. 7. Spectral functions in the (0, 0) — > (0, n) direction 
for points separated by an equal distance. 



FIG. 8. Spectral functions in the (tt, 0) — » (tt, tt) direction 
for points separated by an equal distance. 



FIG. 9. Spectral functions in the direction (0, tt) — > (tt, tt) 
for points separated by an = equal distance. 



FIG. 10. Spectral functions in the direction (0, 0) — > (tt, tt) 
for points separated by an = equal distance. 



FIG. 11. Quantum mechanical probability w that an addi- 
tional hole is in the polaron state at a particular site in the 
stripe structure for a wave function related to a dominant 
peak in the spectral function at the energy about — 1.60t for 
p = 3D(tt,0). 



FIG. 12. Quantum mechanical probability w that an ad- 
ditional hole is in the polaron state at the site n in the 
stripe structure for a wave function related to a recognizable 
peak in the spectral function at the energy about — 2.23t for 
p = 3D(tt, 0). The hole stripes are in the th and 4 th column. 



FIG. 13. Quantum mechanical probability w that an addi- 
tional hole is in the polaron state at a particular site in the 
stripe structure for a wave function related to a dominant 
peak in the spectral function at the energy about — 1.63t for 
p = 3L>(0,tt). 



FIG. 14. Spectral functions in the (0, 0) — * (tt, 0) direction 
for points separated by an equal distance. Calculation was 
performed for the structure depicted in Fig. |H| 



FIG. 15. Spectral functions in the direction (0, 0) — > (tt, tt) 
for points separated by an equal distance. Calculation was 
performed for the structure depicted in Fig. |E| 
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